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Abstract. Intermittency, measured as log (F(r)/3), where F(r) is the flatness of velocity increments at scale 
r, is found to rapidly increase as viscous effects intensify, and eventually saturate at very small scales. This 
feature defines a finite intermediate range of scales between the inertial and dissipation ranges, that we shall 
call near-dissipation range. It is argued that intermittency is multiplied by a universal factor, independent 
of the Reynolds number Re, throughout the near-dissipation range. The (logarithmic) extension of the 
near-dissipation range varies as \/\ogRe. As a consequence, scaling properties of velocity increments in the 
near-dissipation range strongly depend on the Reynolds number. 
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1 Introduction 

Statistics of developed turbulence are commonly investi- 
gated by means of (longitudinal) velocity increments 5v(r) 

" across a distance, or scale, r. At r w Lq, where Lq rep- 
resents the characteristic scale of the stirring forces (the 
integral scale of turbulence), fluid motions are statistically 
independent and the probability density function (pdf) of 
Sv(Lq) is found nearly Gaussian. At smaller scales, in- 
trinsic non-linear fluid dynamics operate and turbulent 
motions become intermittent; fluid activity comes in in- 

■ tense locally-organized motions embedded in a sea of rel- 
atively quiescent and disordered eddies (see |1I2| for first 
numerical indications). As a consequence, the pdf of Sv(r) 
develops long tails and becomes strongly non-Gaussian. 
Deviations from the Gaussian shape may be quantified by 
the flatness, defined as 



For a centered Gaussian distribution F = 3; as long tails 
develop F increases. F(r)/3 may therefore be roughly 
thought of as the ratio of intense to quiescent fluid motions 
at scale r. In that sense, we shall assume in the following 
that log (-F(r)/3) provides a quantitative measure of in- 
termittency |3] . 

The normalized (to the Gaussian value) flatness is plot- 
ted as a function of the scale ratio t/Lq for two turbulent 
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flows in Fig. ^ Experimentally, a particular attention has 
been paid to the size of the hot-wire probe and to the 
signal-to-noise ratio; small-scale velocity fluctuations are 
expected to be suitably resolved 0]. The "instantaneous 
Taylor hypothesis" (see [HJ for details) has been used to 
estimate spatial velocity increments and reduce modula- 
tion effects |3I4| . Details about the (standard) numerical 
integration of the Navier-Stokes equations can be found in 
UJ. We observe at scales r > Lq, F(r) ~ 3, in agreement 
with the picture of disordered fluid motions: There is no 
intermittency, since the flatness F(r) is independent of the 
scale r and (almost) equal to the Gaussian value F = 3. 
At smaller scales, F(r) displays a power-law dependence 
on r: Intermittency grows up linearly with log(l/r). This 
scaling behavior is inherent to the inertial (non-linear) 
fluid dynamics and refers to the so-called inertial range. 
The exponent £f ~ —0.1 is found very consistent with al- 
ready reported values for homogeneous and isotropic tur- 
bulence 0. Interestingly, F(r) exhibits a rapid increase as 
viscous effects intensify, and eventually saturates at very 
small scales. This rapid increase of intermittency, which 
occurs over a range of scales that we shall call the near- 
dissipation range, is the main concern of this article. We 
shall here argue how intermittency in the near-dissipation 
range is related to the build-up of intermittency in the 
inertial range; the Reynolds-number dependence of this 
phenomenon will be also addressed. 

There has been a considerable amount of works on 
intermittency in the inertial range (see |5] for a review). 
Dissipation-range intermittency has received much less at- 
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Fig. 1. Scale dependence of the flatness of longitudinal ve- 
locity increments for two different turbulent flows: (solid line) 
a turbulent jet (R\ — 380) [Hj and (dashed line) a direct nu- 
merical simulation (R\ = 140) 0. The control parameter Re* 
is defined as Re/R,, where Re denotes the (usual) Reynolds 
number; the empirical constant R* ~ 56 (10111) . Appendix El 
is devoted to the (modified) Reynolds number Re,. 



tention. In 1967, Kraichnan conjectured "unlimited inter- 
mittency" for the modulus of velocity Fourier modes at 
very high wavenumbers |12| . Although no proof was ex- 
plicitly established for the Navier-Stokes equations, Frisch 
and Morf provided in 1981 strong mathematical argu- 
ments (occurrence of complex-time singularities) [5] in 
support of Kraichnan's conjecture. Following a previous 
study carried out by Paladin and Vulpiani in 1987 |13| . 
Frisch and Vergassola suggested in 1991 that multifractal 
(local) exponents h of velocity increments, Sv(r) ~ r , 
are successively turned-off as viscous effects intensify |14| : 
As the scale r decreases, only the strongest fluctuations 
(low h) survive while the others are extinguished by the 
viscosity. This mechanism reinforces the contrast between 
intense and quiescent motions, and thus provides a phe- 
nomenological explanation for the increase of intermit- 
tency in the near-dissipation range. However, as the re- 
maining intense motions concentrate on a smaller and 
smaller fraction of the volume, this approach again pre- 
dicts "unlimited intermittency" for the velocity increments 
Sv(r), in the limit of vanishing scale r. 

As mentioned above, our experimental and numerical 
data indicate that intermittency, measured by the flat- 
ness of velocity increments, does exhibit a blow up in the 
beginning of the dissipation range but eventually satu- 
rates in the far-dissipation range. At this point, it should 
be mentioned that a spurious limitation of intermittency 
may stem from a lack of accuracy (or resolution) in ve- 
locity measurements or numerical simulations. However, 
a special care has been taken here to reduce this effect 0]. 
In the following, we will argue that the observed satura- 
tion of intermittency is not an artefact, but refers to some 
peculiar properties of turbulence at very small, dissipative 
scales. 



2 A multiplicative cascade description of 
intermittency 

In the present study, the issue of intermittency in the dissi- 
pation range is reconsidered. The saturation of the flatness 
in the limit of vanishing scale r is recovered by assuming 
that the velocity field is smooth (regular) in quiescent-flow 
regions (as already suggested in jlftllfij ): 5v(r) is not zero 
but behaves as r in these regions; Sv(x, r) ss rd x v(x). This 
key assumption is here recast in a multiplicative approach 
of velocity-increment statistics along scales, as brought 
forward by Castaing et al. in rT7j. We shall then demon- 
strate that it is possible to gain quantitative results, with- 
out ad hoc parameters, on dissipative-range intermittency: 
The amplification of intermittency in the near-dissipation 
range and the extension of the near-dissipation range are 
explicitly estimated as a function of the Reynolds number. 

The build-up of intermittency along the whole range 
of excited scales is related to the distortion of the pdf of 
5v(r). In order to account for this distortion, let us for- 
mally introduce a random independent multiplier (3(r / Lq), 
connecting the statistics of 5v(r) at scales r and Lq: 

Sv{r) = /3(f) x 5v(L ) for r < Lq . (1) 
Lq 

The integral scale L is taken as the reference scale. Eq. ([Ql 
should be understood in the statistical sense, i.e., the pdf 
of 5v(r) equals the pdf of (3(r/L ) x Sv(L ) (see [IB] for a 
Markovian description). The multiplier /3(r/Lo) is consid- 
ered as a positive random variable. This approach there- 
fore restricts to |<5w(f)| or to the symmetric part of the pdf 
of Sv(r); the skewness effects are beyond the scope of the 
present description. 

From Eq. |TjQ), it can be established that 

where P r and GV,l denote respectively the pdf of 6v(r) 
and of \og(3(r/Lo). The pdf of 5v(Lq) may be considered 
as Gaussian, as mentioned in the introduction. Once Pl 
is known, P r is fully determined by G r ,L - 

The so-called propagator kernel G t ,l |4I17I19I2()] is 
characterized by the whole set of coefficients C n (r/Lo), 
defined by the expansion 

/ °° n \ 

<09(t-)) p > = e*P E C n(-r)-r for all p. (2) 



By construction, C n (r/Lo) is the n -order cumulant of 
the random variable log/3(r/Lo). For our purpose, we shall 
only focus on the first two cumulants: The mean 

d(f ) = <log/?(f- )) (3) 

and the variance 

C 2 (f ) ee (log 2 /3(f)) - (log/3(f )) 2 . (4) 

Lq Lq Lq 
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Experimentally, higher-order cumulants of log (3{r / Lo) 
are found very small compared to C\(r/Lo) and C^r/Lo) 
|19I20| . This motivates our main interest in the mean and 
variance of log (3(r/L ) [2]. However, the exact shape of 
the propagator kernel G r ,L is n °t relevant for the follow- 
ing analysis: Our arguments apply to the mean and the 
variance but does not require Gr,L being Gaussian. This 
point should be unambiguous. We shall demonstrate that 
considering the mean and the variance of G Tt L is valuable 
in order to describe the amplification of intermittency in 
the near-dissipation range. 

What can be said about C\{r/L ) and Ci{rjLo)l First 
of all, it is straightforward to get from Eq. Q and Eq. (j2J) : 

( oo n \ 

n=l ^° 

by assuming that 5v(Lq) is a zero- mean gaussian variable 
of variance a 2 (i.e., (\5v(L )\ p ) — K p a p ). 
— In the inertial range, 



G r , Lo (log/3) 



(Nr)r>=v 



(5) 



This is the postulate of universal power-law scalings [9122] . 
It follows that Ci(r/Lo) and C^ijfLo) behave as linear 
functions of log(r/io) : 

T T 

Ci(— ) = cilog( — ) and 

V V 

C2\—) = C2 log( — ) in the inertial range, (6) 
Lo Lo 

where c\ and C2 are universal constants 19j. The depar- 
ture from the Kolmogorov's linear scaling law C, v — p/3 is 
directly related to c 2 < 0: 

P 2 

C P = cip + c 2 — + ■■■ 

In our framework, the build-up of intermittency (along the 
inertial range) is related to the increasing width of G rj £ 
with the decreasing scale r, stating that the second-order 
cumulant C 2 (r/L ) increases as r decreases. 

- In the far-dissipative range, velocity increments are 
proportional to the scale separation r, which leads to 

d(f- ) = log(f ) + Cf ss - and 

Ci[—— ) = C2 1SS ' in the far-dissipation range. 
Lo 

The constants Cf ms - and C^ 3 ' a priori depend on the 
Reynolds number, here defined as 



Re 



(7) 



where L is the integral scale pointed out by Eq. ©, a 
denotes the standard deviation of Sv(L ) and v is the kine- 
matic molecular viscosity. 

- Finally, the inertial-range and far-dissipation-range 
behaviors of Ci(r/L ) and GiirjLo) match in the near- 
dissipation range. We will see that this matching is (very) 
peculiar. 



inertial 




far — dissipation 



logQ3ite.^) 



Fig. 2. Sketch of the distortion of the pdf of log f3(r/Lo) along 
scales r, respectively in the inertial, near-dissipation and far- 
dissipation ranges. The width of the pdf increases rapidly when 
crossing the near-dissipation range from right to left; with de- 
creasing scale r. 



3 The near-dissipation range 

Following Paladin and Vulpiani one considers that 
viscous effects at a given scale r only affect the fluctu- 
ations of 5v(r), for which the local Reynolds number is 
smaller than a certain constant In the classical phe- 
nomenology of turbulence i?* is fixed to unity |9| , but for 
our purpose, i?* is kept as an empirical constant. 

In our multiplicative cascade framework, the previous 
hypothesis writes 



< R* 



This condition is equivalent to 
/3(f )Re* ff 



< 1 



(8) 



where i?e* denotes the (modified) Reynolds number 

R^ = ^ ■ (9) 



The subscript * indicates that i?* is not a priori equal to 
unity; this point is clarified in the Appendix 1X1 

The propagator kernel G r .L is sketched in Fig. [21 for 
various scales r (the integral scale Lo is fixed) as a function 
of log(/3i?e*r/Lo): G r ,L moves from right to left as the 
scale r decreases. At a given scale r, fluctuations j3(r/Lo) 
which satisfy Eq. |(HJ undergo dissipative effects. As viscos- 
ity strongly depletes these affected fluctuations, a signifi- 
cant stretching of the left tail of G r ^ Q is expected (see Fig. 
EJ. Viscous effects then lead to a strong distortion of G r ,z, ; 
a significant increase of Ci(rlLo) follows. This argument 
provides a qualitative explanation of the increase of inter- 
mittency due to (non-uniform) viscous effects. Within this 
picture, the near-dissipation range may be viewed as the 
range of scales r marked by the entering of G r x a m t ne 
viscous domain, and the leaving of G r ,L from the inertial 
domain (see also Fig- EJ ■ 

In order to pursue a more quantitative analysis, an 
explicit definition of the near-dissipation range is required. 
To do so, let us introduce the two characteristic scales rj + 
and rj- (see Fig. given respectively by 



log(Re 



L ' 



C 2 ( 



Lo 



(10) 
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Fig. 3. The near-dissipation range is given by r\- < r < 
77+. The characteristic scales 77 + and r]~ are defined by Eqs. 
HUB and mil : The scale r/ + may be viewed as the smallest 
limiting scale for which the propagator G r ,L is not affected by 
viscous effects; only non-linear dynamics prevail at r > n + ■ The 
scale rf~ may be viewed as the largest limiting scale for which 
the propagator G r ,L does not undergo non-linear effects; only 
viscous dynamics prevail at r < rf~ . 



and Ci(^)+log( J Re*^' 



(id 



According to our previous considerations, 77+ may be seen 
as the scale marking the entering of G r ,L m the viscous 
domain, and rj- as the scale marking the leaving of Gr,i 
from the inertial domain. In other words, r\- < r < rj + 
may be considered as the near-dissipation range; this will 
be our (explicit) definition. 

It is natural to match the inertial-range and dissipation- 
range behaviors of C\{r/L Q ) at the characteristic scale 77, 
for which the propagator kernel, approximately centered 
around its mean value, extends equally over the inertial 
and viscous domains (see Figs. |21 and |3J. This statement 
writes 



= 



and yields 



d(f ) + io g (ife.f ; 



= 



By assuming that the intermittency correction on c\ is 
very small, so that c\ w 1/3 according to Kolmogorov's 
theory (22], one obtains that 77 coincides with the notorious 
Kolmogorov's scale tjk, based on the modified Reynolds 
number Re*: 



T] = 

Furthermore, Ci lss - 
and l|T2*|). it follows 



m = L Q (Re*y 3 ^ . (12) 
= 1/2 log Re*. From Eqs. 0, JTTJ 



ygg) 



|log(at) + io g jte, 
-21og(fe)-§logite, 



21og(^) ' 



and by considering that G„^ results from the distortion 
of G n+ .L , one obtains (see Appendix [5] for a kinematic 
proof) 

2,/c 2 (^+)«|log(^+) (13) 
L 3 ?7_ 



2 i /C 2( |-)«21og ( ^) 



(14) 



These results finally yield log(?7 + /7/) = log(?7/77_) and 



V = Vv+v- ■ 

In logarithmic coordinates the dissipative scale 77, given by 
Eq. I|12(l . lies at the center of the near-dissipation range: 77 
separates the inertial range and the far-dissipation range, 
as originally proposed by Kolmogorov. This is a first result 
concerning the near-dissipation range. 
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Fig. 4. Scaling behavior of C2(r/Lo) for the turbulent jet 
and the numerical simulation. In the inertial range, Ci{r) = 
C2 log(r/Z/o) with C2 ~ —0.025. The two dissipative scales 
77+ and 77_, defined by Eqs. 1)100 and 1)11(1 . are indicated for 
both flows. The "amplification law" between Ci(j)- / Lq) and 
6*2(77+ /-^o) is very well satisfied. We observe also that the Kol- 
mogorov's dissipative scale 77, given by Eq. 1)120 . lies approxima- 
tively at the center of the near-dissipation range: 77- < r < 77+ . 



4 The amplification law 



From the above computation, one can derive the "ampli 
fication law" 

Lq 4 Lq 



(15) 



which characterizes the increase of intermittency in the 
near-dissipation range. The 9/4 factor relies on the (rea- 
sonable) approximation that c\ ~ 1/3. This does not mean 
at all that intermittency is ignored in our approach; it is 
just assumed here that the value of the parameter c\, en- 
tering in the description, can be considered very close to 
its Kolmogorov value. Anyhow, the same reasoning could 
be pursued by keeping c\ as a free parameter. In that case, 
the multiplicative factor would express as 4/(1 + ci) 2 . Ex- 
perimentally, one finds c% = 0.37 ± 0.01 |17| . 
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Interestingly, the amplification of intermittency in the 
near-dissipation range is universal, independent of the 
(very high) Reynolds number. Let us also remark that this 
"amplification law" may serve as a useful benchmark for 
the (experimental or numerical) resolution of the finest 
velocity fluctuations; one should be able to differentiate 
(true) viscous damping and spurious filtering (see |23| for 
such a debate). 



A C 2 (r/L ) 



log i?e* 




log(?7_/Lo) log(?7 + /Lo) 

logfa/Lo) = -3/41ogfle, 



log(r/L ) 



Vlog Re* 




Cx{r/L Q ) 



Fig. 5. Sketch of the scaling behavior of Ci(r / ' Lq) and 
C2(r/Lo) at high Reynolds number. 



5 A unified picture of intermittency 



At very high Reynolds number, yj C^t/Lq) becomes neg- 
ligible compared to log -Re* in the near-dissipation range. 
It follows from Eqs. and the "amplification law" 

that rj + w rj- ~ r\ at leading order in logi?e*. This yields 



point, it should be emphasized that r\+ can not be assim- 
ilated to the Taylor's microscale A. Indeed, log(ry + /?/) oc 
\/log Re* while on the contrary log(A/rj) oc logi?e». 

Velocity increments follow a log-infinitely divisible law 
|24I25| in the inertial range, since all cumulants are pro- 
portional to a same function of the scale — C p (r/Lo) = 
c p log(r / Lq) for all p — but this log-infinitely divisibility 
can not pertain in the near-dissipation range, according 
to the sketch in Fig. |SJ 

Intermittency has been related to log(F(r)/3) in the 
beginning. Using Eq. (TJJ, one can derive the exact equa- 
tion 



log 



F(r) 



(17) 



where C p (r/Lo) is thep-th order cumulant of \og(3(r/Lo). 
At leading order, this yields 



log 



(18) 



which suggests that \og(F(r)/3) and C(r/L ) should be- 
have in a very comparable way. Fig. ^ and Fig. 01 are in- 
deed very similar. One may thus consider that C2(r/L ) 
provides an alternative measure of intermittency, less in- 
tuitive than log (F(r)/3) but physically more tractable (as 
demonstrated by this study). A specific test of Eq. i|18|) is 
provided in Fig. [H] 




0.1 0.2 

log(F(r)/3)/log(Re) 



C 2 (^)-C 2 (^)cxlogi?e, 
and Gl (^)-a(2t)«-|]og(2t) 

L Lq 3 



with log(— ) cx v/log Re* . (16) 
V- 



Fig. 6. We observe that Eq. 1181 is well satisfied in the inertial 
range (for the turbulent jet and the numerical simulation). A 
small departure is observed at small scales, when viscous effects 
intensify. This departure may be attributed to the growing of 
cumulants of order p > 3, neglected in Eq. I|r81 , 



These behaviors are sketched in Fig. |SJ One obtains 
that the (logarithmic) extension of the near-dissipation 
range (~ ^log Re*) becomes negligible compared to the 
extension of the inertial range (~ logi?e„). This is consis- 
tent with the tendency observed in Figs. Q and 0| At this 



Before concluding this study, let us mention that it is 
quite direct to generalize the previous analysis to iV-order 
velocity increments: 

S' l 'v(r) — v(x + r) — v(x), 
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S^v(r) = v(x + r) 



2v(x H — ) + v(x), etc. 



r^ p are preserved 



Inertial-range scalings (\S^ N ^v{r)\ p ) 
(as argued in [2]]) but (\S^v(r)\ p ) ~ r N ? in the far- 
dissipative range. As a result, the definition of the near- 
dissipation is unchanged but the "amplification law" be- 
comes 



^0 



9 / (N 
4 [~ 



1) 



cf)(ft) 



(19) 



The amplification factor depends on the order TV of the 
velocity increment. This feature allows us to discriminate 
the inertial range and the near-dissipation range: Inertial- 
range scalings do not depend on N, while on the opposite, 
near-dissipation-range scalings depend drastically on N. 
Finally, the amplification factor in Eq. 1)19)1 diverges with 
N, tending toward Kraichnan's view of unlimited inter- 
mittency (for the fluctuations of velocity Fourier modes) 
as N -> oo [E]. 



6 Conclusion 

A unified picture of velocity- increment intermittency, from 
the integral scale to the smallest (excited) scales of mo- 
tion, is proposed. It is explicitly stated how far-dissipation 
range and inertial-range intermittencies match in the near- 
dissipation range. Especially, a universal "amplification 
law" determines how intermittency of velocity gradients 
is linked to the build-up of intermittency in the inertial 
range. The results are found in good agreement with our 
experimental and numerical observations. 

Beyond these precise results, this study indicates that 
there are some peculiar and interesting physics around the 
Kolmogorov's dissipative scale of turbulence. Such issue 
may be of great importance, for instance, in the modelling 
of mixing properties of turbulence, which mainly rely on 
the behavior of gradient fields |26|, 

Finally, we would like to insist on the fact that this de- 
scription leads to predictive results which could be used as 
tests for the suitable resolution of (very) small-scale fluc- 
tuations, to distinguish "probe effects" and true viscous 
damping. Relations between this study and the so-called 
property of Extended Self- Similarity [27J should deserve 
some interests as well. 
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A The (modified) Reynolds number Re* 

The empirical constant which is abusively fixed to 
unity in the classical phenomenology of turbulence 9 , 
may be linked to the Kolmogorov's constant ck (see [25] 
and references therein). 



In Kolmogorov's 1941 theory, ck can be defined through 
the second-order velocity structure function: 

{ {5v {r)?)=c K (ef\ 2 '\ 

where (e) denotes the mean dissipation rate. Here, inter- 
mittency corrections are obviously omitted. By the use of 
Eq. one can write 



(Mr)) 



2/3 



which yields 



ck 



(20) 



In this monofractal description, the near-dissipation range 
is degenerate and reduces to the Kolmogorov's scale tjk- 
The second-order moment of velocity gradient expresses 
as 

. 2 . _ ((Sv( VK )) 2 ) 



(vk) 2 



(21) 



By assuming homogeneous and isotropic turbulence, the 
mean dissipation rate writes (e) = lbv((d x v) 2 ) . By com- 
bining the Eqs. 1)20)) and 1)21)1 . together with the definition 
of t)k given by Eq. one gets 



ck 



R„ 



2/3 



— ) or _R* = 15c 



3/2 



(22) 



Eq. 1)221) indicates that i?* is eventually much greater than 
unity. Following |10| . the empirical value R„ rs 56 corre- 
sponds to ck ~ 2.4. This value is in good agreement with 
experimental and numerical estimations ck ~ 2 |28| . 



B Kinematic proof of Eqs. (fT5j ) and (fT4j ) 

We provide a kinematic proof of Eqs. (|13fl and 1)141) with 
the help of Fig. □ 

- on the one hand, one derives from Eq. i|l(J|) that the 
distance 

x A - xa> = ' 2 \J C2 (~]~) ■ 

— on the other hand, the "position" of the point M, de- 
fined by the standard deviation around the mean, moves 
from A to A' within the inertial domain (see Fig. [7J with 
a typical "velocity" 



dx M 



dlogr/Lo 



(ci + 1), 



as r decreases from rj+ to r}-. Within this representation 
(see [23 for details), the variable log r/L may be viewed 
as "time" . The correction due to the change of width of 
G r ^L is neglected. Indeed, this correction expresses as 
\ /2yJ C2 / \og{r / Lq) and can therefore be omitted in the 
near-dissipation range. One then gets 



xa - XA' 



>g(^ 
3 rj. 



) by taking c\ 
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G r , i0 (log/3) 



viscous domain 



inertial domain 
r — rj+ 




log(/3ite*f-) 



Fig. 7. When the scale r decreases, the propagator kernel 
Gr,L moves from right to left (a denotes the standard devia- 
tion). At scale r = r]+, the points A and B are defined by the 
standard deviation around the mean. As r decreases from 77+ 
to r\- , the points A and B move respectively to A' (within the 
inertial domain) and to B' (within the viscous domain). 
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